Upper Limit on the Cosmological Gamma-ray Background 
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We show that the current extragalactic gamma-ray background (EGB) measurement below 100 
GeV sets an upper limit on EGB itself at very high energy (VHE) above 100 GeV. The limit is 
conservative for the electromagnetic cascade emission from VHE EGB interacting with the cosmic 
microwave-to-optical background radiation not to exceed the current EGB measurement. The cas- 
cade component fits the measured VHE EGB spectrum rather well. However, once we add the 
contribution from known source classes, the Fermi VHE EGB observation exceeds or even violates 
the limit, which is approximated as E 2 dN/dE < 4.5 x 10~ 5 (£/100 GeV) -0 ' 7 MeV/cm 2 /s/sr. The 
upper limit above 100 GeV is useful in the future to probe the EGB origin and the new physics like 
axion-like particles and Lorentz-invariance violation. 



I. INTRODUCTION 

The origin of the unresolved extragalactic diffuse 
gamma-ray background (EGB) radiation has been a big 
puzzle in astrophysics and astroparticle physics. The 
EGB was first discovered by the SAS-2 satellite 
EGRET (Energetic Gamma-Ray Experiment Telescope) 
on board the Compton Gamma-Ray Observatory con- 
firmed the EGB spectrum at 0.03-50 GeV Q. Recently, 
LAT (Large Area Telescope) on board the Fermi Gamma- 
ray Space Telescope [Fermi] made a new measurement 
of the EGB spectrum from 0.2 to 100 GeV The ob- 
served integrated EGB flux (E > 100 MeV) is 1.03 x 10~ 5 
photons/cm 2 /s/sr with a photon index of 2.41±0.05. 
This power-law spectrum extends up to 600 GeV based 
on the very recent preliminary EGB spectrum reported 
by the Fermi collaboration Q ■ 

EGB is composed of various unresolved gamma-ray 
sources. Point sources detected by EGRET and Fermi 
are guaranteed to contribute to EGB. Those are namely 
blazars [e.g. radio galaxies^, starburst galaxies [e.g. 
0, H|, high latitude pulsars [J and gamma-ray bursts 
(GRBs) It is expected that blazars, radio galaxies, 
and starburst galaxies explain 22.5 ± 1.8% @, 25±H% 
[| and 4 - 23% of the unresolved EGB, respectively. 
Other extragalactic sources have also been discussed as 
the origin of EGB, although they are still not detected 
in gamma-ray [see |6|, and references therein] . 

Very high energy (VHE; > 30GeV) gamma-rays prop- 
agating through the universe experience absorption by 
the interaction with the extragalactic background light 
(EBL) via electron-positron pair production [e.g.[Tl|. As 
discussed in @ , if the EGB radiation originates from cos- 



mological sources, the EBL absorption signature should 
appear in the spectrum above ~30 GeV. However, the 
measured EGB spectrum shows a single power-law up to 
600 GeV [1, 3] . This may pose a serious problem for the 
current models. 

Electron-positron pairs created by VHE gamma-rays 
with EBL scatter the cosmic microwave background 
(CMB) radiation via the inverse Compton (IC) scatter- 
ing and generate secondary gamma-ray emission compo- 
nent (the so-called cascade emission) in addition to the 
absorbed primary emission [e.g. 12]. At redshift z, the 
scattered photon energy E~ iC appears at lower energy 
than the intrinsic photon energy E^^, typically 
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The cascade component is also expected to contribute to 
EGB QJ, G3]. Recently Murase et al. (2012) [3] con- 
strained the cosmic energy density of gamma-rays using 
the cascade component contribution to the EGB from 
the Fermi measurement. 

In this paper, we generalize the argument in a conser- 
vative way for arbitrary cosmological sources, and set 
an upper limit on EGB by itself with the new Fermi 
EGB data, in particular, on VHE EGB by requiring the 
cascade emission not to exceed the currently observed 
EGB below ^30 GeV. Taking into account the guaran- 
teed source's contributions, we find that the current EGB 
measurement already self-limits the VHE EGB fluxes 
above 100 GeV to PeV as 



odN * ( E 

E dE <4 - 5Xl ° lOOGeV 



-0.7 



'Electronic address: 
t Electronic address: 



yinoue@slac.stanford.edu 



kunihito.ioka@kek.jp 
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for cosmological sources, which is inconsistent with the 
current EGB measurement by Fermi. We further dis- 
cuss the requirements for possible origins of VHE EGB. 
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Hereafter we use the standard cosmology (Ho, Qm, ^a) 
(70.0 km/s/Mpc,0.3,0.7). 



II. THE EGB SPECTRUM 

The EGB spectrum in the 
MeV 2 /cm 2 /s/sr/MeV is calculated as 
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where E ^ s is the observed photon energy, c is the 
light speed, t is the cosmic time, \dt/dz\ = Hq(1 + 
z)y/flM(l + z) 3 + Ha, and T 77 (-E bs, z) is the gamma-ray 
opacity for -E b s from z. We assume z max = 5, which does 
not affect our results. 

The comoving volume emissivity dj/dE 7 (Ej,z) 
[ph/s/MeV/cm 3 ] is given by the intrinsic plus cascade 
emission, j = j int + j cas . 

The intrinsic emission can be characterized by a few 
parameters, 



^(E„z) = h^(l + zY 
dE^ 7 ' ' \ 0, 



E < E 

Ej > -Z?max j 



(4) 



where E 1 is the photon energy in the rest frame, the 
spectral shape is a power-law with a photon index F p h 
and an cutoff at E max , the ^-evolution is given by /3 cvo , 
and jo is the normalization. 

To give a conservative upper limit, we adopt the EBL 
model by as shown in Figure [TJ The EBL intensity 
by ll| is close to the galaxy counts which are the lower 
limit of the EBL. As the EBL density becomes higher, 
the EGB upper limit gets tighter. 



III. CASCADE EMISSIVITY 

Following [HE3, 

we calculate the cascade emissivity 

djcas/dE^ as: 
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where tic(z) is the energy- loss time of an electron with a 
Lorentz factor j e and mass m e by the inverse Compton 
(IC) emission in the local rest frame, 

tic(z) = , 3meC , ; « 7.7X10 13 (*) 1 (1+z)- 4 s, 

(7t is the Thomson scattering cross section, and hcmb(^) 
is the CMB energy density at z. We consider the CMB 
photons only here, since the EBL energy density is two 
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FIG. 1: The EBL model fll] used in this paper is shown 
by solid line. For comparison, other EBL models [l5l | are 
shown as indicated in the figure. The integrated brightness 
of galaxies [l(J (minimum EBL; filled symbols) and current 
measurements of the EBL [TtJ (open symbols) are shown as 
indicated in the figure. The upper limit from TeV gamma-ray 
observation by MAGIC is also shown by solid line with 
arrows. 



orders magnitudes lower than that of CMB. die/d'ye is 
the electron injection spectrum: 
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and d 2 N~ e e / dtdE~ is the scattered photon spectrum per 



unit time by the IC scattering: 



d 2 N^ 
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with f(x) = 2x\n(x) + x + 1 - 2x 2 , (0 < x < 1) and 
x = E^/A-f^e. Here, E lt i — 2-f e m e c 2 is the energy 
of intrinsic photons and dncMB /de is the CMB photon 
density. The integration region over the Lorentz fac- 
tor, 'Yg, is 7e,min ^ 7e ^ 'fe.m&x: l^max — H max /2lTL e C 

and 7e !m in = (E 1 /e) 1 ^ 2 /2. Since the cooling time tic is 
usually shorter than the comoving time, we assume that 
pairs generate photons at the pair creation site. Because 
of this fast cooling, the low energy photon spectrum be- 
low 100 MeV becomes r pri = 1.5. We do not take into 
account this spectral effect since it does not affect the 
VHE spectrum. 

We iteratively calculate Eq. ([5]) by substituting 
(2dj cas /d'ye) (1 — e _T "' 7 ) for dj e /dj e in order to include 
IC scatterings due to pairs from reabsorption of cascade 
photons [Tg| . 

The intergalactic magnetic field (IGMF) effect is not 
important in this study. Although IGMF bends motion 
of created charged pairs and some fraction of beamed 
emission is lost, off- axis sources complement this loss. 
The synchrotron cooling is also not effective for pairs 
typically created outside galaxies. 
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FIG. 2: Upper limit on EGB by requiring the cascade emis- 
sion not to exceed the EGB data below 100 GeV (cascade- 
limit case) in the model-independent way. The observation 
is still consistent with this upper limit. We set the photon 
index r p h = 1.5, the z-evolution index fi cvo — 0.0, and the 
maximum energy E mBX = 60.0 TeV in Eq. @. Dotted, dot- 
dashed, double dot-dashed and dashed curves show the in- 
trinsic spectrum (no absorption), absorbed, cascade, and total 
(absorbed+cascade) EGB spectrum, respectively. Thick solid 
curve with arrows show the upper limit. The filled square 
points show the observed EGB data with the 11-months Fermi 
data [3j. The circle points show the observed EGB data with 
the 24-months Fermi data Error bars represent l-cr un- 
certainty of the data. 



IV. RESULTS 

There are three main parameters: the photon index 
F p h; the z-evolution index f3 cvo ; and the maximum energy 
-Emax in Eq. (@|. As shown later, r p h » 1.5 gives the most 
conservative limit on EGB. r p h = 1.5 is also expected 
to be the hardest photon index in the diffusive shock 
acceleration [201 ] . 

Most of astrophysical sources show positive cosmolog- 
ical evolution (f3 evo > 0) with the comoving number den- 
sity increasing with redshift. High-frequency-peaked BL 
Lacs (HBLs), elliptical galaxies, and cluster of galaxies 
are known to show no or negative evolution, /3 evo = 
(HI, /3 evo = -0.86 [Hj], and /3 evo = -1.0 [23, respec- 
tively. Hereafter, we study (3 cvo < to minimize the 
EBL effect, otherwise we notice it. 



A. Model-Independent Self-Limitation method 

Figure [5] shows a typical example of the self-limitation 
method. Given (r ph , /3 cvo , £ max ) = (1.5,0.0,60 TeV), we 
can only adjust the normalization of EGB by requiring 
the cascade emission not to exceed the data observed by 
Fermi Q (cascade-limit case). Then we obtain the upper 
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FIG. 3: Upper limits on the EGB spectrum for various pho- 
ton index F p h parameters are shown here. Dashed, dotted, 
and dot-dashed curves correspond to r ph = 1.0, 1.5, and 2.0, 
respectively. Thick solid curves represents the upper limit for 
each parameters. For the case of T p h = 2.0, we show the cas- 
cade component by double dot-dashed curve. The upper limit 
on EGB with r p h = 2.0 is set by requiring the primary com- 
ponent not to exceed the observation (primary-limit case). 
The upper limits on EGB for r p h = 1.0 and 2.0 are tighter 
than the case of r p h = 1.5. The observational data are the 
same as those in Fig. [2] 



limit above 100 GeV, well approximated as 

dN f E \ ^ 

E 2 < 1.1 x 10~ 4 [ ) MeV/cm 2 /s/sr. 

dE VlOOGeVy ' 1 ' 

This is still consistent with the current observation. The 
normalization of the EGB upper limit is determined by 
the observed EGB data at ~10 GeV. Although ~60 TeV 
emission from extragalactic sources has not been ob- 
served yet, we adopt E max = 60 TeV here to constrain 
the VHE EGB. We show the cases with different max- 
imum energies and spectral models later. Here we do 
not include the guaranteed source classes' contribution 
to EGB. As shown below, the limit violates the observa- 
tion once we take into account the known source's con- 
tribution such as flat spectrum radio quasars (FSRQs), 
BL Lacs, radio galaxies, and starburst galaxies. 

Figure O is the same as Figure [5] but for T p h = 1.0, 
1.5, and 2.0. Even if we increase the normalization of 
the input EGB or take T p h < 1.5 to explain the VHE 
EGB data, the cascade flux increases at the same time. 
Then the limit becomes stronger than the case with T = 
1.5. The model with r p h = 2.0 shows a typical example 
that EGB is limited by the primary component (primary- 
limit case). Softer input spectrum results in a stronger 
upper limit. This is because the VHE EGB spectrum is 
determined by the absorbed component alone, not by the 
cascade component, although the cascade spectral shape 
is almost independent of the primary spectrum. Thus, 
Tph ~ 1-5 is the most conservative case. 
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FIG. 4: Upper limits on the EGB spectrum for various param- 
eters are shown here. We show the total (absorbed + cascade) 
spectrum with the photon index F p h = 1.5. Dot-dashed, dot- 
ted, and dashed curves correspond to the maximum energy 
£ max = 0.6, 6.0, 60 TeV for /3 cvo = 0.0, respectively. Double 
dot-dashed, and triple dot-dashed curves are the limits for 
/3 CVO =-2.0, and 2.0 for £ max = 60 TeV, respectively. Thick 
solid curves with arrows represent the upper limits for each 
parameter sets. The observational data are the same as those 
in Figure [2] 



FIG. 5: Upper limits on the EGB spectrum for various 
spectral models with /3 cvo = 0.0 are shown here. Dashed, 
dotted, and dot-dashed curves correspond to (r p h,-Emax) = 
(1.5,60 TeV), Mrk 421-like spectrum, and Blackbody spec- 
trum with Ebb = 1 TeV, respectively. Thick and thin 
curve represents the intrinsic spectrum and the total (ab- 
sorbed+cascade) spectrum, respectively. Thick solid curves 
represent the upper limit for each parameter. The observa- 
tional data are the same as those in Fig. [2] 



Figure @] shows the upper limits on EGB for r p h =1.5 
with various /3 evo and E max . Negative evolution with 
/?cvo < eases the upper limit on EGB. In the neg- 
ative evolution case, the dominant EGB contribution 
comes from low redshift and the absorption effects are 
small. Since a large fraction of VHE emitting sources 
locates inside the gamma-ray horizon without suffering 
the EBL attenuation, the contribution of cascade emis- 
sion becomes minor. The maximum energy as low as 
Em ax — 0.6 TeV also eases the upper limit. This is be- 
cause the cascade emission appears only at ~ 0.3 GeV fol- 
lowing Eq. [T] However, there are no known sources that 
have a large contribution to EGB and a spectral cutoff 
at ^TeV. For example, imaging atmospheric Cherenkov 
telescopes (IACTs) detects > 0.6 TeV emission from sev- 
eral nearby blazars [e.g. 24]. 

Each curve in Fig. Ogives the EGB upper limit in each 
energy band for a fixed /3 ovo because the limit is basically 
set by the original flux before absorption at the maxi- 
mum energy. For /3 evo = 0, the limit is approximated by 
Eq. Note that T p h ~ 1.5 is the most conservative 
case as discussed above. 

With /3 ovo = 2.0, i.e. positive evolution, the upper 
limit violates the EGB measurement for (r p h,-E m ax) = 
(1.5, 60 TeV). To avoid the inconsistency with the mea- 
sured EGB, the VHE emissivity beyond the gamma-ray 
horizon should be turned off or less than that inside the 
horizon. Therefore, if the EGB origin is cosmological, the 
source might have hard spectra and show a no or nega- 
tive cosmological evolution. Recent Fermi analysis shows 
that FSRQs have (3 = 5.7 (25[. No or negative-evolution 



sources reported in gamma-ray or in other wavelength are 
HBLs 21], elliptical galaxies [22] and clusters of galaxies 
[23l ]. Gamma-ray emission from latter two has not been 
confirmed yet and is not likely enough for EGB at least 
in the case of cluster of galaxies [26[ . 

Fig. [5] shows the upper limits for non-power-law spec- 
tral models in the case of no evolution. Here we show a 
Mrk 421-like spectrum and a blackbody shape spectrum. 
For a Mrk 421-like spectrum, we use a log-parabola func- 
tion as 



■(E~/, z) oc 



Ebi 



-rph+siogc-; 



x (l + z)' 3 - (10) 



where we use the best fit parameters for Mrk 421 as Ebr = 
0.3 TeV, r ph = 2.48, and 5 = 0.33 [H- For a blackbody 
shape spectrum, we adopt 



■(Ej, z) cx 



exp( J B 7 /£ , BB) - 1 



x (1 + zf° 



(11) 



where we set -Ebb = 1 TeV. The upper limits for both 
of a Mrk 421-like spectrum and a blackbody spectrum 
models comes lower than that for T = 1.5 spectral model. 
Therefore, T « 1.5 is the most conservative case even if 
we consider these non-power-law spectral models. 

Fig. [6] shows the upper limit for (r p h,/3 e vo, E max ) = 
(1.5,0.0,60 TeV). Here we show the contribution from 
each redshift ranges. At VHE region, only sources at 
z < 0.5 can contribute to the EGB due to the EBL sup- 
pression. This means that the dominant VHE EGB con- 
tribution at each energy roughly comes from inside the 
gamma-ray horizon by EBL attenuation. Therefore, we 
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to explain the EGB below 10 GeV by guaranteed sources 
as in Fig. the EGB measurement violates the upper 
limit above 100 GeV. The limit is approximated by Eq. 
[5] If we take /3 cvo = —4, the upper limit becomes con- 
sistent with the measured spectrum with 1-sigma differ- 
ence. However, there are no known sources showing such 
a strongly negative evolution in any wavelengths. In the 
case of E max = 6 TeV, f3 cvo < —6 is required. Therefore, 
the VHE emissivity at redshift z > 0.5 should be low. 
On the other hand, when we try to make the upper limit 
consistent with the VHE EGB data as in Fig. [SJ the total 
EGB contribution from cosmological sources is ^2 sigma 
below the measured EGB below 10 GeV. In this case, the 
limit is approximated by 1.6 times Eq. [2J 



FIG. 6: Upper limits on the EGB spectrum for various red- 
shift bins are shown. Thick solid curve with arrows shows the 
upper limit. Solid, dotted, dot-dashed, double dot-dashed, 
triple dot-dashed and dashed curves show < z < 0.01, 
0.01 < z < 0.1, 0.1 < z < 0.5, 0.5 < z < 2, 2 < z < 5, 
and < 2 < 5 EGB spectrum, respectively. 



can constrain the VHE emissivity of the universe at each 
energy by future VHE EGB measurements. 



B. Self-Limitation method with known sources' 
contribution 

There are guaranteed source classes that contribute 
to EGB detected by EGRET or Fermi. It is expected 
that blazars, radio galaxies, and starburst galaxies ex- 
plain 22.5 ± 1.8% 0, 25±?|% @ and 4 - 23% Q of the 
unresolved EGB, respectively. Then, - 70% of EGB will 
be explained by known source classes. We need to sub- 
tract them to evaluate the VHE EGB upper limit, since 
the residual is the only room for the cascade plus ab- 
sorbed emission from VHE EGB. 

For FSRQs, we adopt the model by [25| (hereafter 
MA12). For BL Lacs, we use the model by [5| (hereafter 
FermilO). Since the EBL absorption effect is not taken 
into account in [5j , we include the EBL attenuation model 
11| by assuming f3 evo = for BL Lacs. For radio galaxies, 
we use the model by @ (hereafter YI11). For starburst 
galaxies, we use the IR luminosity function model from 
Q (hereafter SV11) and the power law model from Q 
(hereafter Fermil2). We renormalize the SV11 model by 
a factor of 0.8 to avoid the total (FSRQ+BL Lacs+radio 
galaxy+starburst galaxy) contribution exceeding the ob- 
served EGB data, or the VHE upper limit becomes zero. 

Figs. [JJ and [3] shows the upper limit on the EGB tak- 
ing into account the known sources' contributions where 
we use the SV11 model and the Fermil2 model for star- 
burst galaxies, respectively. Here we show the case of 
(r ph ,/3 ovo ,£; max ) = (1.5,0.0,60 TeV). The upper limit 
on EGB is derived from the sum of the VHE EGB cas- 
cade and guaranteed sources' contribution. When we try 



V. DISCUSSION AND CONCLUSION 

There are a few possible scenarios to explain VHE 
EGB. (i) Sources have hard spectra T p h ~ 1.5 with a cut- 
off at 60 TeV and a strongly negative evolution f3 evo < — 4 
in Eq. (fj|. If a cutoff is at 6 TeV, /3 cvo < —6 is required, 
(ii) More transparent EBL generates weaker absorption 
effect and eases the limit. However we use the EBL 
model close to the minimum (integrated flux of galax- 
ies). Even if we use it, the upper limit is still below 
the observation, (iii) Pair production process could be 
affected by new physics such as Lorentz-invariance viola- 
tion (2?| and axion-like particles [28j]. (iv) Dark matter 
annihilation/decay in the local group [2^, [3(| can avoid 
the EBL absorption effect, (v) Sources only contribut- 
ing to EGB at < 10 GeV may complement the resid- 
ual between data and model in Fig. [5] High latitude 
pulsars and radio-quiet AGNs are possible candidates. 
First, pulsars observed in gamma-ray have a cutoff at 
~ 5 GeV [3l[ . Second, radio-quiet AGNs may contribute 
to the EGB at < 10 GeV (H, although Fermi does not 

34j |. If non-thermal electrons 



see radio-quiet AGNs H, 



exist in a corona above the accretion disk, a power-law 
tail will appear in hard X-ray and gamma-ray band [32| . 
(vi) The EGB measurement has uncertainties. The EGB 
is deduced by subtracting the foreground emission from 
our Galaxy which is still not fully understood. For ex- 
ample, in the analysis in p|, the Fermi bubble [35j is not 
subtracted. 

In the scenario (i), we can reject known Fermi gamma- 
ray source classes as the origin of VHE EGB. First, 
blazars and radio galaxies detected by Fermi do not show 
negative cosmic evolution 0, [|| . Secon d, g amma-ray ob- 
served galaxies do not show r p h = 1.5 [3y|. Even if they 
can create such hard spectra, TeV emission is internally 
absorbed by the interstellar radiation [37j • Third, pul- 
sars will not contribute at VHE band as discussed above. 
Although a power-law tail in the VHE band have been 
recently reported for the Crab pulsar [38| , the photon in- 
dex r p h = 3.8 is softer than that of the observed VHE 
EGB r ph = 2.41. 

One of the most likely source classes is TeV selected 
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FIG. 7: Upper limit on the EGB spectrum by adding 
the known sources' contribution for (r p h, Pevo, E max ) — 
(1.5,0.0,60 TeV). The dotted, dot-dashed, double dot- 
dashed, and triple dot-dashed curves correspond to contri- 
bution from FSRQ [H, MA12], BL Lac [|, FermilO], radio 
galaxy 0, YI11], and starburst galaxy @, SV11], respectively. 
For BL Lac, we include EBL absorption effect [rH ] by as- 
suming /3 cvo — 0. For starburst galaxy, we use the SV11 IR 
luminosity function model renormalized by a factor of 0.8. 
The thin solid curve show the cascade+absorbed contribu- 
tion from VHE EGB (self-limit method). The dashed curve 
show the total emission of FSRQs, BL Lacs, radio galaxy, 
starburst galaxy and our model. The thick solid curve with 
arrows shows an upper limit which is approximated by Eq. [2] 
The observational data are the same as those in Figured] 



HBL (TeV HBL) which is not detected by Fermi but 
by the current IACTs (39[. Although their cosmological 
evolution is still unknown, their spectrum is hard and 
could have a cut-off at ~ 10 TeV. Interestingly these hard 
TeV emission can be explained by high energy cosmic-ray 
induced intergalactic cascade [40|. A new IACT array 
Cherenkov Telescope Array (CTA) [4l| is expected to 
detect > 100 blazars including TeV HBLs g2|. CTA 
will enable us to statistically study their evolution and 
contribution to VHE EGB. 

Low luminosity (LL) GRBs may also explain VHE 
EGB, although LL GRBs have not been detected in 
gamma-ray. LL GRBs might show a negative cosmo- 
logical evolution since LL GRBs are only discovered 
at low redshift [e.g. [43[. The total energy budget of 
LL GRBs is large enough to explain Ultra-High En- 
ergy Cosmic Rays (UHECRs) 44]. The UHECR in- 
tensity is ~ 1.0 x 10" 5 MeV/cm 2 /s/sr at ~ 10 19 eV 
45] which is comparable to the observed EGB intensity 
~ 3.0 x 10" 5 MeV/cm 2 /s/sr at 600 GeV Q. 

It is important to measure VHE EGB more precisely 
and at higher energy, such as by Fermi, CTA and CALET 
(CALorimetric Electron Telescope) [4fJ. Our upper limit 
above 600GeV will help future measurements to unveil 
the EGB origin. Just by detecting the EGB above TeV, 
we can put a meaningful lower limit on the number den- 



FIG. 8: Same as Fig. but we use the Milky Way model 
by @, Fermil2] for starburst galaxies. The thick solid curve 
with arrows shows an upper limit which is approximated by 
1.6 times Eq. d 

sity of the EGB sources because a source should reside 
in the gamma-ray horizon that is small at high energy. 
For example, if CTA measures the EGB at 60TeV, the 
gamma-ray horizon is ~ 40 Mpc and the EGB source 
number in the entire sky should be larger than N m - m = 
4?ri ; EGB(60TeV)/FcTA(60TeV)_- 6[F EGB (60TeV)/3 x 
10~ 7 MeV/cm 2 /s/sr] [see alsoflil where we assume 
CTA sensitivity as F CT A(60TeV) = 1.0 x 10- 12 erg/cm 2 /s 
[4~fl ] . We note that the total emissivity within the gamma- 
ray horizon may have a large dispersion (such as Poisson 
fluctuation) at high energy because of the small horizon 
size, which may lead to a violation of our upper limits. In 
other words, the number of sources within the gamma- 
ray horizon may be larger at high energy region than our 
expectation due to the local distribution fluctuation. 

Anisotropy in the VHE EGB above 30 GeV is an alter- 
native key to understand VHE EGB, since the monopole 
peak depends on sources [48|. The anisotropy below 50 
GeV has already been investigated [49j . 

The effect of IGMF to the cascade component is not 
critical in this study. However, with weak strength of 
IGMF, Fermi may have already detected the cascade 
emission alone from some sources and classified them into 
unassociated sources [50j . If their contribution had been 
already subtracted from the EGB, the upper limit on 
the VHE EGB would be eased. If the cascade emission 
alone had been detected, the absorbed VHE flux would 
be brighter than the Fermi's vF v sensitivity. Then, the 
current IACTs should have detected those VHE emission, 
since their sensitivity in vF v in the VHE region is compa- 
rable to that of Fermi in the GeV region. However, TeV 
HBLs have not been detected by Fermi yet [3!| . There- 
fore, there would be no sources whose cascade component 
alone is resolved. 

Strong IGMF would also ease the upper limit on the 
VHE EGB. When the IGMF strength is above 3.26 /xG, 
the magnetic energy density takes over the CMB en- 
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ergy density in the Thomson regime. Then, the syn- 
chrotron cooling is more effective than the IC cooling. 
The strong IGMF will suppress the cascade emissivity in 
gamma-ray. It is, however, also known that the mag- 
netic field strength of lobes of AGN is typically 1/iG 
[Ell ]. Thus, the strength of IGMF may be smaller than 
1//G with the scale of AGNs' lobes - lOOkpc dJ. This 
scale is shorter than the mean free path of pair creation 
- 20 Mpc(n E BL(z)/0.1cm- 3 ). Therefore, strong IGMF 
would not affect the cascade emissivity in gamma-ray sig- 
nificantly 

The upper limit is also applied to the EGB produced 
by the UHECRs via the intergalactic cascade [see e.g. 
[13]. In order to explain the VHE EGB by the UHECR 
cascade with known sources' contribution to EGB, a neg- 
ative evolution may be necessary for the UHECR sources. 

In this paper, we develop a new method to constrain 
the cosmological EGB by using the EGB itself. VHE pho- 
tons propagating the universe are absorbed by EBL and 
create electron-positron pairs. Created pairs generate 
secondary gamma-ray emission via the IC scattering of 
the CMB photons. We constrain the VHE EGB by com- 
paring this regenerated emission with the current EGB 
measurement. Our method also provides upper limits 
on EGB above 600GeV for future observations, such as 
Fermi, CTA and CALET. Our self-limits are useful for 
identifying the origin of EGB. 



We show that the current EGB measure- 
ment sets an upper limit of E 2 dN/dE < 
4.5 x 1CT 5 (£;/100 GeV)~ - 7 MeV/cm 2 /s/sr to the 
VHE EGB above 100 GeV from cosmological sources, 
where we take into account the known sources' contri- 
butions. The current EGB measurement by Fermi [ij 
violates the predicted upper limit. In order to make con- 
sistent with the observed EGB data in the cosmological 
origin scenario, possible origins should show strongly 
negative cosmological evolution, hard spectrum, and a 
cut-off at ~ 10 TeV. This kind of sources, however, has 
never been reported yet neither in gamma-ray nor in 
other wavelength. 
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